Coherent combining pulse bursts in time domain

ABSTRACT

A beam combining and pulse stacking technique is provided that enhances laser pulse energy by coherent stacking pulse bursts (i.e. non-periodic pulsed signals) in time domain. This energy enhancement is achieved by using various configurations of Fabry-Perot, Gires-Tournois and other types of resonant cavities, so that a multiple-pulse burst incident at either a single input or multiple inputs of the system produces an output with a solitary pulse, which contains the summed energy of the incident multiple pulses from all beams. This disclosure provides a substantial improvement over conventional coherent-combining methods in that it achieves very high pulse energies using a relatively small number of combined laser systems, thus providing with orders of magnitude reduction in system size, complexity, and cost compared to current combining approaches.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 61/918,485, filed on Dec. 19, 2013. The entire disclosure of the above application is incorporated herein by reference.

GOVERNMENT CLAUSE

This invention was made with government support under Grant No. DE-SC0006377 awarded by the U.S. Department of Energy, and N6601-11-1-4208 awarded by the U.S. Defense Advanced Research Projects Agency. The government has certain rights in this invention.

FIELD

The present disclosure relates to lasers and more specifically to coherent combining of multiple laser beams and temporal pulse stacking to increase pulse energy.

BACKGROUND

In general, beam combining of multiple lasers allows overcoming power and energy limitations of each individual laser. Beam combining of multiple laser signals can be currently achieved using a variety of methods such as active coherent phasing, incoherent spectral combining, passive self-locked combining, or incoherent spatial addition of multiple laser beams. In each of these methods, each parallel laser channel generates an identical signal and all parallel output signals are then combined with the total combined power proportional to the sum of individual powers from all N laser channels. The maximum achievable power can never exceed this total sum, as fundamentally limited by the power conservation law. These beam combining methods can be applied to continuous wave signals as well as pulsed signals. When these currently used methods are used to combine pulsed beams, the combined energy per each pulse can never exceed the sum of individual pulse energies from all the channels.

In case of pulsed signals, it is possible to achieve simultaneous pulse combining and time-domain energy redistribution such that in a combined beam pulse repetition rate is reduced and, therefore, combined energy per pulse now increases proportionally to both this repetition rate decrease and to the total sum of individual pulses from all the channels. This increases combined pulse energy faster than a linear proportionality to a number of parallel channels N. As described in applicant's previous work, this increase could be proportional to N² of the number of channels, thus significantly reducing combined laser array size and complexity for reaching the same combined pulse energy as with the above-described methods. Further details regarding this previous work can be found, for example in U.S. patent application Ser. No. 14/403,038 which is incorporated in its entirety herein. The technique described in this previous work is applicable only to periodically pulsed signals. Redistribution of energy between the pulses in time domain requires that a beam-combining element would provide a time delay longer than the time duration over which pulse energy is redistributed. This time delay is always associated with the size of the combining element. For example, in the N² combining method described in applicant's previous work, a Fabry-Perot or any other configuration of a resonant cavity can be used, with its round-trip length L=c·ΔT, where ΔT is the each-channel periodic-signal repetition period, as shown in FIG. 1. This means that pulse repetition rate at the input of the system has to be high, in order to be compatible with a practically short length of the combining cavity. This means that achieving high energy at high repetition rates is associated with impractically high average powers. For example, a 3 m long roundtrip cavity corresponds to 100 MHz repetition rate of the periodic pulse input. Assuming that each pulse in the input periodic signal is at about 1 mJ of energy (expected energy of ˜1 ns pulse at the maximum of each parallel amplification channel), then the average power per each amplified beam from each parallel channel is approximately 100 kW.

This section provides background information related to the present disclosure which is not necessarily prior art.

SUMMARY

This section provides a general summary of the disclosure, and is not a comprehensive disclosure of its full scope or all of its features.

In one aspect, a method is provided for coherently combining pulses of a pulsed optical signal in the time domain. The method includes: receiving an optical input signal comprised of a periodic pulse train; splitting the input signal into two optical signals; forming at least one cluster of pulses in each of the two optical signals, where amplitude and phase of pulses in each cluster varies amongst pulses; and coherently combining the two optical signals in the time domain using a resonant cavity and thereby forming an optical output signal having a solitary pulse for each cluster of pulses, where the two optical signals are combined in a manner such that pulse energy of the solidary pulse exceeds N times highest energy of any given pulse in the cluster of pulses from the two optical signals.

Clusters of pulses may be formed by amplitude modulating and phase modulating the two optical signals. In particular, the clusters of pulses in the two optical signals are formed complementary with each other in the spectral domain. The two optical signals may be combined using a Fabry-Perot interferometer.

This method may be extended to coherently combine pulses from N optical signals. The method includes: receiving an optical input signal comprised of a periodic pulse train; splitting the input signal into N optical signals; amplitude modulating each signal in the N optical signals to form clusters of pulses in each signal, such that amplitude of pulses in each cluster varies amongst pulses; phase modulating each signal in the N optical signals, such that phase of pulses in each cluster varies amongst pulses; and coherently combining each of the phase modulated signals in the time domain using one or more resonant cavities and thereby forming an optical output signal having a solitary pulse for each cluster of pulses, where the N optical signals are combined in a manner such that pulse energy of the solidary pulse exceeds N times highest energy of any given pulse in the cluster of pulses from the N optical signals.

Prior to coherently combining the phase modulated signals, each of the phase modulated signals may be amplified. Additionally, the N optical signal can be phase locked.

In some embodiments, the phase modulated signals are coherently combined using a plurality of Fabry-Perot interferometer cavities arranged in either a parallel cascading manner or a sequential cascading manner.

In another aspect of this disclosure, a method is provided for coherently combining (or stacking) pulses from a single optical signal in the time domain. The method includes: receiving an optical input signal comprised of a periodic pulse train; forming at least one cluster of pulses in the pulse train, where amplitude and phase of pulses in the at least one cluster varies amongst pulses; and coherently combining pulses in the at least one cluster in the time domain to form an optical signal with a solitary pulse using a resonant cavity, where the resonant cavity has an entirely reflective cavity. The resonant cavity may be further defined as a Gires-Tournois interferometer cavity.

Pulse energy can be further increased through the use of N parallel channel amplification array. This method includes: receiving an optical input signal comprised of a periodic pulse train; amplitude modulating the optical input signal to form at least one clusters of pulses, such that amplitude of pulses in each cluster varies amongst pulses; phase modulating the optical input signal, such that phase of pulses in each cluster varies amongst pulses; splitting the optical input signal into N optical signals; amplifying each of the N optical signals; coherently combining each of the amplified optical signals into a combined signal using a beam combiner; and coherently combining pulses in the at least one cluster of pulses in the combined signal in the time domain using a resonant cavity having an entirely reflective cavity, thereby forming an optical output signal having a solitary pulse for each cluster of pulses.

Pulses from the pulse clusters can be coherently combined using a Gires-Tournois interferometer. In some embodiments, N Gires-Tournois interferometers can be arranged in a sequential cascading manner, such that the energy increase of a solitary pulse is approximately 2N. In other embodiments, Gires-Tournois interferometers can be arranged in two or more sequential stages, where each stage includes N Gires-Tournois interferometers arranged in a sequential cascading manner and the Gires-Tournois interferometer in each subsequent cascade are of equal roundtrip length, which is at least 2N times longer than in the previous cascade.

Further areas of applicability will become apparent from the description provided herein. The description and specific examples in this summary are intended for purposes of illustration only and are not intended to limit the scope of the present disclosure.

DRAWINGS

The drawings described herein are for illustrative purposes only of selected embodiments and not all possible implementations, and are not intended to limit the scope of the present disclosure.

FIG. 1 is a diagram illustrating a conventional beam combining technique;

FIG. 2 is a diagram illustrating a coherent beam combining and pulse stacking technique;

FIG. 3 is a diagram depicting an example embodiment of a beam combining and pulse stacking system;

FIGS. 4A-4C are diagrams depicting example implementations for resonant cavities used for beam combining and pulse stacking;

FIGS. 5A and 5B are diagrams depicting example arrangements for cascading resonant cavities;

FIG. 6A illustrates a specific example for sequentially cascading an array of Fabry-Perot interferometer cavities;

FIGS. 6B-6G are graphs pulse profiles for six beams input into the cascading arrangement of FIG. 6A;

FIG. 7 is a diagram illustrating that the pulse stacking technique can be extended to a resonant cavity with only one input mirror;

FIG. 8 is a diagram depicting an example embodiment of a pulse stacking system using fully reflective resonant cavities;

FIG. 9 is a diagram depicting an alternative embodiment of a pulse stacking system using fully reflective resonant cavities;

FIGS. 10A-10C are diagrams depicting example implementations of fully reflective resonant cavities;

FIGS. 11A and 11B are diagrams illustrating pulse stacker arrangements which may be used in the embodiments of FIGS. 8 and 9;

FIG. 11C-11G are graphs depicting examples of calculated pulse-bursts which can serve as input to the sequential cascade of FIG. 11B;

FIG. 12 is a diagram depicting an arrangement for a pulse-burst stacker having two cascaded stages;

FIG. 13 is a diagram depicting another arrangement for a pulse-burst stacker having two cascaded stages;

FIG. 14 is a diagram depicting coherent pulse stacking with a Gires-Tournois interferometer;

FIG. 15 is a diagram depicting amplitudes of all incident, in-cavity and output pulses in a conceptual 2-mirror GTI cavity;

FIG. 16A is a diagram illustrating m cascaded equal-roundtrip GTI cavities;

FIG. 16B is a diagram illustrating m×m multiplexed different-roundtrip GTI cavities;

FIGS. 17A-17C are graphs depicting calculated intensities and phases of input pulse bursts and output pulse intensities of four cascaded cavities; and

FIGS. 18A-18C are graphs depicting calculated intensities and phases of input pulse bursts and output pulse intensities of four-by-four multiplexed cavities.

Corresponding reference numerals indicate corresponding parts throughout the several views of the drawings.

DETAILED DESCRIPTION

Example embodiments will now be described more fully with reference to the accompanying drawings.

Instead of dealing with periodic signals, this disclosure presents beam combining and pulse stacking (in time domain) of solitary pulse bursts. This approach allows decoupling of the resonant cavity length from the pulse repetition rate, since this cavity length now determines pulse repetition rate within the solitary pulse burst, while repetition rate between these solitary pulse bursts is completely unconstrained by the combiner size, and, therefore, can be selected completely arbitrarily.

FIG. 2 illustrates the beam combining and pulse stacking technique using a resonant cavity 20 with two input beams, one beam 22 configured to be reflected and the other beam 23 transmitted by the cavity. The resonant cavity 20 is characterized by the fact that the two end mirrors defining it's cavity are each partially reflecting (i.e. each of their reflectivity is less than 100%). In an example embodiment, the resonant cavity 20 is further defined as a Fabry-Perot Interferometer cavity. Other types of resonant cavities are also contemplated by this disclosure.

The concept is best illustrated by considering this resonant cavity combiner 20 in reverse. That is, if a solitary pulse 24 is incident into the resonant cavity 20, it will split into two beams, one corresponding to a reflection-port input 26, another to a transmission-port input 27. This “reversed” pulse now will produce a time-decaying burst of pulses in each of the two beams. In general, each pulse in each of the bursts is characterized by pulses of varying amplitude and phase, which can be easily calculated. Pulse repetition rate in each of the burst will be identical and equal to the RC combiner round-trip time. Consequently, if the “reversed” sequence of each of the pulse bursts is produced for the input reflection and transmission channels (which includes reversing the order of the pulses and replicating correct amplitudes and phases for each of the pulse in the sequence), then a single beam is produced at the output of the RC combiner with all the pulses stacked in time domain, such that a solitary pulse containing the total energy obtained by summing-up all individual pulse energies from each of the two incident pulse bursts. It can be shown that this pulse energy enhancement with respect to each individual pulse in each of the two input pulse clusters will be larger than the number of channels.

FIG. 3 depicts an example embodiment for implementing the proposed beam combining and pulsing stacking technique which can be contrasted with previous methods. In previous methods, a periodic pulse train is split into N equal-amplitude signals, which are then launched into multiple parallel amplification channels. Each signal in each individual channel are then uniquely phase modulated so that the frequency-comb structure of each periodic pulse train in each channel is shifted in frequency and the phase being different for each channel. Using a suitable selection of these frequency shifts as well as a suitable beam combiner, it possible to reconstitute this periodic pulse train, albeit at an N-times reduced repetition frequency.

Likewise, a periodic pulse train 31 serves as the input to the beam combining and pulse stacking system 30. The periodic pulse train 30 is first split by a beam splitter 32 into N optical signals. In each of the N optical signals, the periodic pulse train is then manipulated to form optical signals having one or more clusters of pulses (also referred to herein as pulse bursts), where amplitude and phase of pulses in each cluster varies amongst pulses. In this example embodiment, amplitude and phase modulation are used to carve-out clusters of pulses from the period pulse train. An optical amplitude modulator 33 and an optical phase modulator 34 are disposed in each parallel channel. The role of the amplitude and phase modulation sequence is to produce pulse burst with the required amplitudes and phases of each individual pulses, so that they can be combined and time-stacked in the properly designed beam combining and pulse stacking arrangements. Amplitude and phase modulation pattern for each of the N signals will in general be different between different signals. It is noted that these clusters of pulses can be repeated in time after some arbitrarily chosen time duration T (which is, however, an integer of the original pulse repetition rate), which is not related, and, therefore, can be much longer, than the combining resonator round-trip time.

Pulses from each of the clusters are combined coherently in the time domain using one or more resonant cavities 37 as will be further described below. It is important that all of the parallel channels be phase locked. Thus, additional phase modulators for phasing-locking the parallel channels (and controlled by properly configured phase-locking circuits) 35 are preferably interposed between phase modulators 34 and the resonant cavity 37. Note that in some configurations, phase modulators 34 and 35 could be combined into one, performing both pulse-burst phase modulation and parallel-channel phasing with the same device. Prior to coherently combining pulses, it may also be beneficial to amplify each of the N optical signals. In this case, an optical amplifier 36 is disposed in each channel between the phase locking circuit 35 and the resonant cavities 37.

Additionally, the resonant cavities 37 should be actively controlled so that each cavity maintains its round-trip time precisely at its prescribed level, with a prescribed phase shift. Note that the RC round-trip time should match the repetition period of the initial pulse train, and, consequently, the pulse repetition period in each pulse burst. However, producing the prescribed phase shift (round-trip phase) for each cavity requires very small (on the order of one optical-cycle of the optical signal) mismatch for each cavity round-trip time. At the output of the N-parallel amplifier array 36, all of the beams should be properly arranged to be incident into the resonant cavities as described in more detail below.

Example implementations for resonant cavity beam combiners and pulse stackers are shown in FIG. 4A-4C. A simple flat-mirror Fabry-Perot interferometer cavity is shown in FIG. 4A. One advantage of this implementation is that the finesse of the combining cavity is generally low, approximately equal to the number of channels combined in a cascaded scheme. Low-finesse cavities are much easier to align than high-finesse cavities. For clarity, “1” indicates input 1 beam, which is configured to be transmitted through the Fabry-Perot interferometer cavity, “2” indicates input 2 beam, configured to be reflected from the Fabry-Perot interferometer cavity; and “3” indicates the combined beam comprised of the transmitted and reflected beams.

In order to avoid significant beam expansion due to propagation distance in a resonant-cavity combiner, it may be beneficial to include curved refocusing mirrors into an Fabry-Perot interferometer cavity as seen in FIG. 4B. The curved mirrors are aligned in such a way that beam size is reproduced exactly every round-trip time in the cavity. Configurations in FIGS. 4A and 4B require external separation between the incident reflection-path and the combined beams. This can be done, for example, using a standard polarization beam splitter and a suitable waveplate or a Faraday rotator. The configuration shown in FIG. 4C has the advantage that all the incident beams are separate from each other and, most importantly, from the combined beam.

To combine more than two optical beams, the resonant cavity beam combiners and pulse stackers need to be cascaded together. Two example arrangements are shown in FIGS. 5A and 5B. In FIG. 5A, Fabry-Perot interferometer cavities may be arranged in a parallel cascading manner. That is, optical beams are paired together (e.g., #1 and #2) and each pair of optical beams are combined together concurrently using a series of cavities. Combined optical beams are likewise paired together (e.g., 51 and 52) and combined with each other using another series of cavities. This arrangement is repeated until two optical beams are combined to form a single optical beam 53.

Fabry-Perot interferometer cavities may also be arranged in a sequential cascading manner as seen in FIG. 5B. In this arrangement, two optical beams (e.g., #1 and #2) are first combined together using one cavity to a form a first combined beam. The combined beam is then combined with a third beam (i.e., #3) using another cavity. This arrangement is again repeated until all of the N optical beams have been combined together to form a single optical beam.

FIG. 6A illustrate a specific example of combining optical beams using a sequential cascading arrangement. In this example, six input beams are combined using an array of five Fabry-Perot interferometer cavities. End-of-cavity mirrors are identical for each cavity. Each cavity is further characterized by a different mirror reflectivity for each of the cavity mirrors, as well as by individually different phase-shifts produced by each resonant cavity. Specific values used in this example are as follows: R ₁=0.7143, Δφ₁=0 R ₂=0.6619, Δφ₂=π/3 R ₃=0.5845, Δφ₃=2(π/3) R ₄=0.4565, Δφ₄=3(π/3) R ₅=0.1956, Δφ₅=4(π/3) It is readily understood that values will differ for other arrangements.

FIGS. 6B-6G depict pulse profiles for the six input beams. More specifically, amplitudes (i.e., peak powers) and phases for each input signal are indicated by red bars and blue circles, respectively. Horizontal axis is labeled in terms of the pulse number in the sequence. Note that earlier time is on the right-hand and later time is toward the left-hand. Specific pulse-burst period value is not important in the calculation and can be freely chosen. For reference, integration of the total pulse energy in each pulse burst (cluster) (red solid line), as well as a horizontal line indicating ⅙ of the total energy obtained by cumulatively summing energies of all the pulses in all the inputs are also shown in the figures. Input each of the six input beams into the sequential cascading arrangement described in FIG. 6A yields an output beam with a single solitary pulse. The important aspect of this example is that the highest (normalized) pulse peak-power amplitude in any input beam/pulse cluster does not exceed 0.08, which means that the total combined-pulse energy is 1/0.08=12.5 times larger than the highest energy pulse in any of the input signals—energy enhancement of more than twice the number of channels.

The beam combining and pulse stacking technique can also be extended to a resonant cavity with only one input mirror as seen in FIG. 7. The resonant cavity can be defined by an input mirror 71 having partial reflectivity and the other mirror 72 defining the cavity having 100% reflectivity. For example, such cavity could be implemented by a Gires-Tournois interferometer. In operation, there is only a single incident beam and a single reflected output beam since all of the power is completely reflected. Thus, no beam combining can be achieved with this cavity. Pulse stacking in time, however, can be achieved similar to the partially-reflecting combiner discussed above. Again, it is possible to calculate the amplitudes and phases for each pulse in the incident pulse burst needed to form a solitary pulse.

FIG. 8 depicts an example embodiment for implementing the pulsing stacking technique using one or more resonant cavities with 100% reflectivity. In this example, a single periodic pulse train 81 is input into a sequence of an amplitude modulator 83 and a phase modulator 84, which again modulate each-pulse amplitude and phase in the solitary pulse burst 85. The pulse-burst can be optionally amplified by an optical amplifier 86 to boost its power, but practical advantages can also be achieved without this amplification. Next, the pulse-burst is incident into a pulse stacker arrangement 87 to form an output beam having a solitary pulse. The pulse stacker arrangement 87 can be implemented by a single fully reflective resonant cavity or multiple fully reflective resonant cavities as further described below.

An alternative embodiment for implementing the pulse stacking technique using one or more fully reflective resonant cavities is shown in FIG. 9. A single periodic pulse train 91 is again input into a sequence of an amplitude modulator 93 and a phase modulator 94, which modulate each-pulse amplitude and phase in the solitary pulse burst 95. In this example, much higher pulse energies and power can be achieved by inserting N parallel-channel amplification array. The solitary pulse-burst 95 is split by a beam splitter 92 into N identical optical signals and thereby forming N parallel channels. In each channel, the optical signal is passed through a phase-locking modulator 96 (each controlled by a properly arranged phase-locking electronic circuit) and an optical amplifier 97. The N amplified optical beams are then coherently combined into a single beam using, for example, conventional phasing schemes and conventional beam combiners as indicated at 98. Example beam combiners may include but are not limited to binary-tree combiners or diffractive optical element (DOE) combiners. Lastly, the recombined and amplified pulse-burst is incident onto a pulse stacker arrangement 99 to form an output beam having a solitary pulse. This embodiment enables achieving pulse energies and powers well beyond the single-channel power limitations. It is also important to note that this N-channel array can be configured for a spectral-coherent combining as described in by Wei-zung Chang, Tong Zhou, Leo A. Siiman, and Almantas Galvanauskas in “Femtosecond pulse spectral synthesis in coherently-spectrally combined multi-channel fiber chirped pulse amplifiers,” Opt. Express 21, 3897-3910 (2013). In this case, this array would allow overcoming not only power limitations, but could also overcome spectral bandwidth limitations of each individual amplification channel. This, for example, can be important when generating ultrashort pulses. In this implementation, the N-channel array should contain spectral beam splitting and combining elements as described in the publication noted above. This array could also be a “hybrid” setup containing a combination of both coherent and spectral-coherent combining arrangements.

FIGS. 10A-10C illustrate example configurations for a fully reflective resonant cavity. In these configurations, the second cavity mirror 101 of the cavity is 100% reflective but otherwise these configurations are analogous to those in FIGS. 4A-4C. Thus, there is only one incident beam marked as “1” and one output beam marked as “2”. A Gires-Tournois interferometer cavity is an example of such fully reflective resonant cavity although other types of resonant cavities are also contemplated.

There are two general constrains on the incident pulse burst, and on the parameters of the resonant cavity. The first constraint is that amplitudes and phases imprinted (using amplitude and phase modulators) on all the pulses in the incident pulse burst should be such that a single stacked pulse is produced at the combiner output. This constraint is relevant both for Fabry-Perot and Gires-Toirnois based combining schemes. This constraint can be formulated mathematically. As an example let's consider combining of two signals with a single Febry-Perot combiner. Let's denote a single stacked pulse at the combiner output as {tilde over (p)}_(s)(t) in the time domain. Spectrum {tilde over (P)}_(s)(v) of this pulse can be obtained by fourier transform of this time-domain signal using a standard fourier-transform definition:

${{\overset{\sim}{P}}_{s}(v)} = {{{??}\left( {{\overset{\sim}{p}}_{s}(t)} \right)} = {\int_{- \infty}^{+ \infty}{{{{\overset{\sim}{p}}_{s}(t)} \cdot {\mathbb{e}}^{{- {\mathbb{i}}}\mspace{11mu} 2\;\pi\;{vt}}}\ {{\mathbb{d}t}.}}}}$ Fabry-Perot spectral transfer functions for reflection {tilde over (F)}^(R)(v) and transmission {tilde over (F)}^(T)(v) are:

${{\overset{\sim}{F}}^{R}(v)} = \frac{\sqrt{R} \cdot \left( {1 - {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \right)}{1 - {R \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}}}$ and ${{\overset{\sim}{F}}^{T}(v)} = {\frac{{\mathbb{e}}^{{\mathbb{i}}\;{\delta/2}} \cdot \left( {1 - R} \right)}{1 - {R \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}}}.}$ In this example we take a Fabry-Perot configuration from FIG. 4 with the two partially-reflecting mirrors equal to each other R₁=R₂=R. Cavity round-trip phase is denoted as δ. We can now express this constraint on the incident pulse burst to produce a single output pulse mathematically as: {tilde over (P)} _(s)(v)·{tilde over (F)} ^(R)(v)+{tilde over (P)} _(s)(v)·{tilde over (F)} ^(T)(v)={tilde over (P)} _(s)(v).

Here {tilde over (P)}_(s)(v)·{tilde over (F)}^(R)(v) describes the spectrum of the pulse burst incident into the reflection port 1 and {tilde over (P)}_(s)(v)·{tilde over (F)}^(T)(v) the spectrum of the pulse burst incident into the transmission port 2, that produce the solitary pulse at the output of the Fabry-Perot combiner. The case with two different reflectivities can be described by using appropriate reflection and transmission transfer functions for such a cavity. Similarly, for the Gires-Tournois pulse stacker case this constraint can be expressed as {tilde over (P)} _(s)(v)·{tilde over (F)} ^(GTI)(v)={tilde over (P)} _(s)(v), where left-hand term defines the incident pulse spectrum required to produce a single solitary output pulse. Here {tilde over (F)}^(GTI)(v) is a Gires-Tournois cavity spectral transfer function:

${{\overset{\sim}{F}}^{GTI}(v)} = \frac{\sqrt{R} - {\mathbb{e}}^{{\mathbb{i}}\;\delta}}{1 - {\sqrt{R} \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}}}$

with R being front-mirror reflectivity.

It is straightforward to generalize this example from a single cavity to multiple cascaded cavities. For cascaded Gires-Tournois cavities the spectral transfer function of the cascade is simply a product of all individual-cavity transfer functions. For cascading Fabry-Perot cavities one needs to use a product of the transmission and reflection transfer functions encountered by each corresponding input-port signal on its path to the output. For example, in the example case in FIG. 9 for the incident pulse burst into the port i=2 this transfer function will be ²{tilde over (F)}^(T)(v)·¹{tilde over (F)}^(R)(v), where ²{tilde over (F)}^(T)(v) is the cavity 2 transfer, and ¹{tilde over (F)}^(R)(v) is the cavity 1 reflection functions.

The second constraint is such that there would be an energy benefit when combining these pulsed bursts. This means that the total energy of the solitary output pulse should be more than N times larger (N being the number of inputs into the combining arrangement) than the highest energy of any individual pulse in all the incident pulse bursts. Otherwise, if this combined-output energy is only N times larger, then this could be achieved with any conventional combining approach. Also, this constraint is relevant only for the Fabry-Perot combining. For the Gires-Tournois combining there always is an energy benefit, as long as a solitary output pulse is produced.

Theoretical description of coherent pulse stacking (i.e. temporal pulse combining) with a single Gires-Tournois Interferometer, from which the requirements for a stacking cavity and input pulse burst parameters can be calculated, is given here.

A reflecting interferometer can be configured either as a linear or a traveling-wave cavity. A linear reflecting cavity is essentially a Fabry-Perot interferometer with one completely reflecting mirror, which is commonly referred to as a Gires-Tournois interferometer (GTI). The practical advantage of a traveling-wave reflecting cavity, shown in FIG. 14, is that it allows one to spatially separate the incident input and reflected output beams. For usage convenience we will also refer to this traveling-wave cavity as a GTI. Let's consider a traveling-wave GTI cavity, consisting of a partially reflecting front-mirror M (with power reflectivity R=r²<1), and K completely reflecting beam-folding mirrors M₁, M₂, . . . , M_(K), schematically shown in FIG. 14. Ideally, we should have R_(k)=1 for all k=1, 2, . . . K, but in practice it will always be R_(k)≈1. We can denote the round-trip cavity transmission as α=r₁·r₂· . . . r_(K), where r_(k) is a corresponding k-th mirror amplitude reflection coefficient. Then coefficient α²=R₁·R₂· . . . . R_(K) describes power loss per each round trip due to the finite reflectivity of the folding mirrors. If the round trip distance in this traveling-wave cavity is P, then the round trip time is ΔT=P/c (here cis the velocity of light), and the round-trip phase is δ=2πP/λ₀=ω₀·ΔT (here λ₀ and ω₀ are, respectively, the signal central wavelength and angular frequency). We can describe electric field transmission through this cavity by a transmission matrix [7]:

$\begin{matrix} {\lbrack T\rbrack = {\begin{bmatrix} 1 & 0 \\ 0 & {\alpha \cdot {\mathbb{e}}^{\mathbb{i}\delta}} \end{bmatrix}.}} & (1) \end{matrix}$

The incident and reflected fields at both sides of the front mirror M can be described by a unitary scattering matrix [S], which can be written in a symmetric form:

$\begin{matrix} {\lbrack S\rbrack = {\begin{bmatrix} r & {it} \\ {it} & r \end{bmatrix}.}} & (2) \end{matrix}$

Here r is the front mirror amplitude reflection and t its amplitude transmission coefficient, which for a lossless (and reciprocal) dielectric mirror are related by t²=1−r², with both r and t being real quantities.

Let's assume that all pulses at the input, output and circulating inside of a GTI cavity have identical complex envelopes {tilde over (p)}_(s)(t), where complex notation accounts for the fact that stacking can be achieved with bandwidth-limited (described by real envelopes {tilde over (p)}(t)≡p_(s)(t)) as well as with chirped (described by complex envelopes {tilde over (p)}_(s)(t)) pulses. As illustrated in FIG. 14, we seek to convert an input pulse burst ^(in){tilde over (p)}(t) into a single output “stacked” pulse:

$\begin{matrix} {{{\,^{i\; n}\overset{\sim}{p}}(t)} = {{{\sum\limits_{j}{{{}_{}^{i\; n}\left. p \right.\sim_{}^{}}(t)}}->{{\,^{out}\overset{\sim}{p}}(t)}} = {{{\overset{\sim}{p}}_{s}(t)} \cdot {{\mathbb{e}}^{{\mathbb{i}}{({{\omega_{0}t} + \phi_{0}})}}.}}}} & (3) \end{matrix}$ Let's choose the time axis reference such that this output stacked pulse is centered at t=0. Then we can enumerate all the pulses in the incident train as

$\begin{matrix} {{{{\,^{i\; n}\overset{\sim}{p}}(t)} = {{\sum\limits_{n = 0}^{\infty}{{{}_{}^{i\; n}\left. p \right.\sim_{}^{}}(t)}} = {\sum\limits_{n = 0}^{\infty}{{{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \cdot {{\overset{\sim}{p}}_{S}\left( {t + {{n \cdot \Delta}\; T}} \right)} \cdot {\mathbb{e}}^{{\mathbb{i}\omega}_{0}{({t + {{n \cdot \Delta}\; T}})}}}}}},} & (4) \end{matrix}$

so that n=0 corresponds to the pulse at t=0 and n=1, 2, . . . are at correspondingly increasing separations n·ΔT at negative t values, as shown in FIG. 14. Here ^(in)Ã_(n)=^(in)A_(n)·e^(iφ) ^(n) is the complex amplitude characterizing the amplitude and phase of the n-th individual pulse ^(in){tilde over (p)}_(n)(t) of the semi-infinite input pulse burst. The pulse repetition period ΔT should be equal to the cavity round trip time ΔT, since interference at the front mirror M occurs only when the incident pulse and the circulating in-cavity pulse arrive at that mirror simultaneously. As shown in FIG. 15, let's denote the amplitude of the circulating pulse arriving at the front mirror M from inside the cavity simultaneously with the ^(in)Ã_(n) input pulse as ^(cav.)Ã_(n). Interference at the front mirror M between these two simultaneously incident pulses ^(in)Ã_(n) and ^(cav.)Ã_(n) produces a “reflected” output pulse with an amplitude ^(out)Ã_(n) and a “transmitted” in-cavity pulse, which after passing mirrors M₁ through M_(K) arrives at M as a subsequent circulating pulse ^(cav.)Ã_(n-1). Amplitudes of all incident, in-cavity circulating, and output pulses are shown in FIG. 15 for a conceptualized 2-mirror GTI cavity. Using the front-mirror scattering matrix [S] given by the Eq. (2), and the cavity transmission matrix [7] given by the Eq. (1) we can describe this interference (plus a subsequent round trip) using the following matrix equation:

$\begin{matrix} {{\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{{cav}.}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {\begin{bmatrix} 1 & 0 \\ 0 & {\alpha \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \end{bmatrix} \cdot \begin{bmatrix} r & {it} \\ {it} & r \end{bmatrix} \cdot \begin{bmatrix} {{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{{cav}.}\left. A \right.\sim_{}^{}} \end{bmatrix}}},} & (5) \end{matrix}$ for each n=0, 1, 2, . . . .

Ideally, when coherent pulse stacking is achieved then ^(out)Ã₀=1 and ^(out)Ã_(n)=0 for all n>0. This means that for all incident pulses prior to the last pulse (i.e. n>0), totally destructive interference between the incident and circulating pulses should eliminate all reflections from the front mirror, thus storing all incident pulses as a single circulating pulse inside the cavity. For the last incident pulse (n=0) totally constructive interference in the reflection direction of that mirror should combine the incident and circulating pulses into a single output pulse, thus extracting all stored circulating energy. Additionally, it is ^(cav.)Ã_(n)=0 and ^(out)Ã_(n)=0 for all n<0, since ideally there should be no field left in the cavity after the last n=0 pulse has passed. Using Eq. (5) above, we can express all these conditions mathematically as:

$\begin{matrix} {{\begin{bmatrix} 1 \\ 0 \end{bmatrix} = {\begin{bmatrix} 1 & 0 \\ 0 & {\alpha \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \end{bmatrix} \cdot \begin{bmatrix} r & {it} \\ {it} & r \end{bmatrix} \cdot \begin{bmatrix} {{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{{cav}.}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = 0},{and}} & (6) \\ {{\begin{bmatrix} 0 \\ {{}_{}^{{cav}.}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {\begin{bmatrix} 1 & 0 \\ 0 & {\alpha \cdot {\mathbb{e}}^{\mathbb{i}\delta}} \end{bmatrix} \cdot \begin{bmatrix} r & {it} \\ {it} & r \end{bmatrix} \cdot \begin{bmatrix} {{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{{cav}.}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = 1},2,{\ldots\mspace{14mu}{\infty.}}} & (7) \end{matrix}$

By solving the linear-equation sets Eq. (6) and, iteratively, Eq. (7), we get all the complex amplitudes of the pulses in the semi-infinite train at the cavity input, and of the corresponding circulating pulses inside the cavity:

$\begin{matrix} \left\{ \begin{matrix} {{{}_{}^{i\; n}\left. A \right.\sim_{}^{}} = r} & \; \\ {{{{}_{}^{i\; n}\left. A \right.\sim_{}^{}} = {{- \left( {1 - r} \right)^{2}} \cdot r^{n - 1} \cdot \frac{{\mathbb{e}}^{{- {\mathbb{i}}}\; n\;\delta}}{\alpha^{n}}}},} & {{{{for}\mspace{14mu} n} = 1},2,3,{\ldots\mspace{14mu}\infty}} \\ {{{{}_{}^{{cav}.}\left. A \right.\sim_{}^{}} = {{- i}\;{t\left( \frac{r}{\alpha \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \right)}^{n}}},} & {{{{for}\mspace{14mu} n} = 0},1,2,{\ldots\mspace{14mu}\infty}} \end{matrix} \right. & (8) \end{matrix}$ Correspondingly, peak power coefficients ^(in)B_(n)=|^(in)Ã_(n)|² of the incident pulses are:

$\begin{matrix} \left\{ \begin{matrix} {{{}_{}^{i\; n}{}_{}^{}} = R} & \; \\ {{{{}_{}^{i\; n}{}_{}^{}} = \frac{\left( {1 - R} \right)^{2}R^{n - 1}}{\left( \alpha^{2} \right)^{n}}},} & {{{{for}\mspace{14mu} n} = 1},2,3,{\ldots\mspace{14mu}\infty}} \end{matrix} \right. & (9) \end{matrix}$ Note that by defining the output stacked pulse using Eq. (3) we chose its peak power coefficient to be normalized to 1: ^(out)B₀=1. In an ideal case, when amplitudes of all the pulses in the semi-infinite input train fulfill the Eq. (8), then all other output pulses are absent, i.e. ^(out)B_(n)=0 for all n≠0. According to Eq. (8), if the GTI cavity round trip phase is chosen to be δ=2πm (where m is an integer) then the last pulse in the input sequence is out of phase with respect to the rest of the pulses in the sequence. It is easy to recognize from Eq. (9) that the peak powers of all the pulses in the input sequence prior to the last pulse (i.e. n=1, 2, . . . ) are described by a decreasing geometrical progression.

Eq. (8) and (9) describe an ideal semi-infinite input pulse sequence, when the only output is the stacked pulse. In practice one needs to truncate this semi-infinite pulse train into a finite pulse burst consisting of N pulses. Since we count pulses from n=0, the very first pulse in this sequence corresponds to n=N−1. Coherent pulse stacking conditions defined by Eq. (6) and (7) still apply to all the pulses in the finite sequence, except for this very first pulse. Indeed, since there are no prior pulses this n=N−1 pulse can only reflect from the front mirror, because it cannot interfere with any prior field in the GTI cavity. Consequently, we can rewrite the amplitude-coefficient condition set by Eq. (7) for this first n=N−1 pulse as

$\begin{matrix} {\begin{bmatrix} {{}_{}^{}\left. a \right.\sim_{N - 1}^{}} \\ {{}_{}^{{cav}.}\left. A \right.\sim_{N - 2}^{}} \end{bmatrix} = {\begin{bmatrix} 1 & 0 \\ 0 & {\alpha \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \end{bmatrix} \cdot \begin{bmatrix} r & {it} \\ {it} & r \end{bmatrix} \cdot {\begin{bmatrix} {{}_{}^{i\; n}\left. a \right.\sim_{N - 1}^{}} \\ 0 \end{bmatrix}.}}} & (10) \end{matrix}$ Here we use the lowercase letters ã instead of Ã to distinguish between the (N−1)^(st) pulse amplitude coefficients in the finite and semi-infinite sequences respectively, as shown in FIG. 15. Consequently, Eq. (8) and (10) give us:

$\begin{matrix} \left\{ \begin{matrix} {{{}_{}^{i\; n}\left. a \right.\sim_{N - 1}^{}} = \frac{{}_{}^{i\; n}\left. A \right.\sim_{N - 1}^{}}{1 - R}} \\ {{{}_{}^{}\left. a \right.\sim_{N - 1}^{}} = {- \left( \frac{r}{\alpha \cdot {\mathbb{e}}^{{\mathbb{i}}\;\delta}} \right)^{N - 1}}} \end{matrix} \right. & (11) \end{matrix}$ From this it is clear that the amplitude coefficient of the first input pulse in the finite sequence is by a factor 1/(1−R) larger than the amplitude coefficient of the same input pulse in the semi-infinite sequence, all the other amplitudes of pulses with n<N−1 are the same in both finite and semi-infinite sequences. The GTI output with a finite sequence input contains both the stacked pulse as well as a weak first-pulse reflection with amplitude ^(out)ã_(N-1). From Eq. (11) describing ^(out)ã_(N-1) it is also clear that this reflection can be made negligibly small by increasing the length N of the incident pulse burst. For example, for a 9-pulse stacking sequence this reflection peak power can be smaller than 10⁻³ of the stacked pulse peak power.

This coherent stacking of multiple pulses into one output pulse containing all the energy of the input pulse sequence is beneficial when amplifying high energy pulses in, e.g. a fiber amplifier, since it enables the amplification of pulses with lower peak powers, thus reducing the detrimental nonlinear effects in an amplifier. This benefit is proportional to the peak-power enhancement factor η=^(out) B ₀/^(max){^(in) B _(n)}=1/^(max){^(in) B _(n)},  (12) where ^(max){^(in)B_(n)} denotes the highest peak-power coefficient in the incident pulse sequence. One can maximize this peak-power enhancement factor by selecting GTI parameters, which minimize the magnitude of ^(max){^(in)B_(n)}. By inspecting Eq. (9) it is straightforward to recognize that ^(max)η is achieved when ^(in)B₀=^(in)B₁. This condition defines the optimum front-mirror reflectivity R_(opt) of a GTI cavity:

${{{}_{}^{i\; n}{}_{}^{}} = {{{{}_{}^{i\; n}{}_{}^{}}->R_{opt}} = {{\frac{\left( {1 - R_{opt}} \right)^{2}}{\alpha^{2}}->{R_{opt}^{2} - {\left( {2 + \alpha^{2}} \right)R_{opt}} + 1}} = 0}}},$ The solution of this quadratic algebraic equation, corresponding to a physically meaningful power reflectivity in the range 0≦R≦1, is

$\begin{matrix} {R_{opt} = {\frac{\left( {2 + \alpha^{2}} \right) - \sqrt{\left( {2 + \alpha^{2}} \right)^{2} - 4}}{2}.}} & (13) \end{matrix}$ Corresponding ^(max)η is

$\begin{matrix} {{\,^{{ma}\; x}\eta} = {\frac{1}{R_{opt}} = {\frac{\alpha^{2}}{\left( {1 - R_{opt}} \right)^{2\;}}.}}} & (14) \end{matrix}$ For an ideal case when α²=1 (all folding mirrors are perfectly reflecting) we have R_(opt)=0.382, and the corresponding highest possible peak-power enhancement for a single GTI cavity of ^(max)η=2.62.

FIGS. 11A and 11B schematically illustrate pulse stacker arrangements 87 which may be used, for example in the embodiments of FIGS. 8 and 9. In FIG. 11A, a single resonant cavity with 100% reflectivity is used stack pulses in the time domain. In FIG. 11B, five resonant cavities are arranged in a sequential cascade. FIGS. 11C-11G depict examples of calculated pulse-bursts which can serve as input to the sequential cascade of FIG. 11B. It is noted that all pulses in the burst are nearly equal in amplitude, and the number of pulses is approximately equal to the double of number of stages in the cascade. This means that energy enhancement X is approximately equal to 2N, where N is the number of stages in the cascade. It is important to note that these equal amplitudes have been achieved by optimizing reflectivities of the front mirrors and each individual stage of the cascade, as well as their individual resonant-cavity phases. Without this optimization, equal amplitudes would not be achieved.

Next, this optimization is described theoretically, showing how stacking of large numbers of equal-amplitude pulses can be achieved by using properly configured sequences of multiple GTI cavities. From the description given earlier, it is clear that when stacking with a single GTI cavity the last “switching” pulse should have its energy comparable to the energy stored in the cavity. However, by arranging multiple GTI cavities in a sequence (or cascade), as shown in FIG. 16A, it is possible to overcome this limitation, and to achieve stacking with an input-pulse burst containing multiple equal-amplitude pulses. We can show this by extending the earlier presented theoretical analysis for a single GTI. Let's use the same notation for cavity parameters and field amplitudes, but add subscripts indicating the corresponding cavity in a sequence, as shown in FIG. 16A. Then Eq. (5) can be rewritten to describe the interference in the k^(th) cavity, occurring between the input, output and circulating fields for the n^(th) input pulse as following:

$\begin{matrix} {\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {{{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{k}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{k}}} \end{bmatrix}\begin{bmatrix} r_{k} & {it}_{k} \\ {it}_{k} & r_{k} \end{bmatrix}}\begin{bmatrix} {{}_{k - 1}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}.}} & (20) \end{matrix}$

Note that now the output field from each cavity becomes the input field for the next cavity in the sequence _(k-1) ^(out)Ã_(n)=_(k) ^(in)Ã_(n). This leads to the following set of equations, which are analogous to Eq. (6) and Eq. (7) from above for the semi-infinite input pulse train case when there are m cavities:

$\begin{matrix} {{{\begin{bmatrix} 1 \\ 0 \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{m}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{m}}} \end{bmatrix}\begin{bmatrix} r_{m} & {it}_{m} \\ {it}_{m} & r_{m} \end{bmatrix}}\begin{bmatrix} {{}_{m - 1}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = {0\mspace{14mu}{and}}}}{k = m}} & (21) \\ {{{\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ 0 \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{k}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{k}}} \end{bmatrix}\begin{bmatrix} r_{k} & {it}_{k} \\ {it}_{k} & r_{k} \end{bmatrix}}\begin{bmatrix} {{}_{k - 1}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = {0\mspace{14mu}{and}}}}{{k = 2},3,{{\ldots\mspace{14mu} m} - 1}}} & (22) \\ {{{\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ 0 \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{1}{\mathbb{e}}^{{\mathbb{i}}\;{\delta 1}}} \end{bmatrix}\begin{bmatrix} r_{1} & {it}_{1} \\ {it}_{1} & r_{1} \end{bmatrix}}\begin{bmatrix} {{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = {0\mspace{14mu}{and}}}}{k = 1}} & (23) \\ {{\begin{bmatrix} 0 \\ {{}_{}^{}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{m}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{m}}} \end{bmatrix}\begin{bmatrix} r_{m} & {it}_{m} \\ {it}_{m} & r_{m} \end{bmatrix}}\begin{bmatrix} {{}_{m - 1}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = 1},2,{{\ldots\mspace{14mu}\infty\mspace{14mu}{and}\mspace{14mu} k} = m}} & (24) \\ {{\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{k}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{k}}} \end{bmatrix}\begin{bmatrix} r_{k} & {it}_{k} \\ {it}_{k} & r_{k} \end{bmatrix}}\begin{bmatrix} {{}_{k - 1}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = 1},2,{{\ldots\mspace{14mu}\infty\mspace{14mu}{and}\mspace{14mu} k} = 2},3,{{\ldots\mspace{14mu} m} - 1}} & (25) \\ {{\begin{bmatrix} {{}_{}^{}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{n - 1}^{}} \end{bmatrix} = {{\begin{bmatrix} 1 & 0 \\ 0 & {\alpha_{1}{\mathbb{e}}^{{\mathbb{i}}\;\delta_{1}}} \end{bmatrix}\begin{bmatrix} r_{1} & {it}_{1} \\ {it}_{1} & r_{1} \end{bmatrix}}\begin{bmatrix} {{}_{}^{i\; n}\left. A \right.\sim_{}^{}} \\ {{}_{}^{}\left. A \right.\sim_{}^{}} \end{bmatrix}}},{{{for}\mspace{14mu} n} = 1},2,{{\ldots\mspace{14mu}\infty\mspace{14mu}{and}\mspace{14mu} k} = 1}} & (26) \end{matrix}$

For each incident pulse n≧0 there are 2m algebraic equations, which completely define all 2m fields _(k) ^(in)Ã_(n) and _(k) ^(cav)Ã_(n) (k=1, m) present for that incident pulse in all m cavities, each as a function of front-mirror reflectivities r₁, r₂, . . . , r_(m)r₁, r₂, . . . , r_(m), round-trip phases δ₁, δ₂, . . . , δ_(m) of all GTI cavities in the sequence, and the pulse number n. We are interested only in the input-field amplitudes into the first GTI cavity ₁ ^(in)Ã_(n), each defined by 2m independent parameters r₁, r₂, . . . , r_(m) and δ₁, δ₂, . . . , δ_(m). Consequently, we can choose to write 2m−1 equations |₁ ^(in)Ã_(n)|²=|₁ ^(in)Ã_(n-1)|² for n=0, 1, . . . , 2m−1, which set the condition for all peak-intensities of the last 2m pulses in the incident pulse burst to be equal among each other. This means that pulse burst with 2m equal-amplitude pulses can be stacked in a sequence of GTI cavities with m reflectivities r₁, r₂, . . . , r_(m), and m−1 round-trip phases δ₁, δ₂, . . . , δ_(m-1) defined by solving these 2m−1 equations. One of the round-trip phases (for example δ_(m)) can be freely selected, and only affects the required individual-pulse phases of the input stacking-pulse burst.

Calculated example for the case of 4 cascaded cavities and the semi-infinite pulse train is presented in FIGS. 17A-17C, showing stacking sequence consisting of eight equal-amplitude pulses with a numerical stacking factor close to 9, since the n=8 pulse in the burst has its peak intensity nearly equal to those of the last eight pulses. The calculated GTI-cascade parameters R_(k)=|r_(k)|² and δ_(k) are as follows: R₁=0.535, R₂=0.526, R₃=0.618, R₄=0.666, δ₁=4.66, δ₂=3.15, δ₃=5.46, R₁=0.535, R₂=0.526, R₃=0.618, R₄=0.666, δ₁=4.66, δ₂=3.15, δ₃=5.46, and we chose δ₄=0. Note that in reality this semi-infinite stacking-pulse train can be truncated to a finite pulse burst, just as in the single-GTI configuration.

Cascaded equal-roundtrip GTI cavities provide multiple-pulse stacking proportional to the number of GTI cavities. It is possible to achieve a substantially quicker increase in stacking factors by using a multiplexed GTI cavity configuration shown in FIG. 16B. This configuration essentially multiplexes a cascade of m equal-roundtrip GTI cavities with an identical cascade of m GTI cavities, but with the second cascade round-trip being more than 2m times longer than that of the first cascade. Numerically calculated example for the case of 4×4 multiplexed cavities and the finite input-pulse consisting of 81 equal-amplitude pulses is shown in FIGS. 18A-18C using the same GTI front mirror reflectivities and cavity round-trip phases as in the previous example of 4 cascaded cavities. It shows that using multiplexing of eight GTI cavities one can achieve stacked peak-power enhancement close to 80.

This multiplexed cascaded scheme can be generalized to the scheme shown in FIG. 12, where the multiplexed sequential cascades contain different number of GTI elements, not just the equal number as described above. The scheme employs two or more stages, each stage comprised of a sequentially cascaded resonant cavities. In this example, a first stage is comprised of N₁ resonant cavities and a second state is comprised of N₂ resonant cavities. These two different stages are distinguished from each other by the fact that the round trip time for all RC cavities in the first cascade is t₁, whereas, the round trip time for all RC cavities in the second cascade is 2N₁ t₁. The input signal now consists of 2N₂ bursts, each burst consisting of 2N₁ pulses. The 2N₁ pulses in each burst are configured in such a way that each burst is stacked into a single pulse at the output of the first cascade. Consequently, at the output of the first cascade, a sequence of 2N₂ pulses is produced. The sequence is configured in such a way (by choosing the suitable phase of each individual pulse-burst in the initial sequence) that the two solitary pulses are in turn stacked again in the second cascade to form an output with a single pulse. Assuming that the equal amplitude pulse format as shown in FIG. 11 is used here the total energy in the stacked pulse is E_(stacked)=2N₂×2N₁×E_(individual). Essentially this represents a quadratic increase with a total number of combining cavities in the cascade. The two stage scheme is merely exemplary as there can be more than two cascades or stages in the sequence.

Signal combining (both spatially and temporally) partially-reflecting resonator cavities and completely-reflecting resonator cavities (temporally combining) can be combined as shown in FIG. 13. Such a configuration represents a general beam combining and pulse stacking arrangement shown in FIG. 3. Note that this is merely an example and different variations of the configuration could be used as well. For example, the parallel cascade can replace the sequential cascade of Fabry-Perot type cavities.

An additional advantageous mode of implementing these techniques is as follows. This pertains to the configuration consisting of a periodic pulse source, from which the required pulse burst is “carved out” using an amplitude and a phase modulators, of an amplifier, and of the Gires-Tournois based pulse combiner. In the prior art there is a technique known as divided pulse amplification (DPA) where solitary pulses from a signal source are first split into a number of pulse replicas using a sequence of optical delay lines, amplified, and then reconstituted back into a solitary pulse using an identical setup of optical delay lines. There is a significant practical advantage of replacing optical delay line pulse splitter in this DPA scheme with the periodic pulse source followed by an amplitude and a phase modulators, as described in this disclosure. With a suitable modulation of amplitudes and phases in the generated pulse burst it could reproduce the effect of the optical delay-line spatial arrangement. But this could be implemented in a much more compact, monolithic in the practice, arrangement, compared to a cumbersome arrangement of spatial delay lines.

The foregoing description of the embodiments has been provided for purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosure. Individual elements or features of a particular embodiment are generally not limited to that particular embodiment, but, where applicable, are interchangeable and can be used in a selected embodiment, even if not specifically shown or described. The same may also be varied in many ways. Such variations are not to be regarded as a departure from the disclosure, and all such modifications are intended to be included within the scope of the disclosure. 

What is claimed is:
 1. A method for coherently combining pulses of a pulsed optical signal in the time domain, comprising: receiving an optical input signal comprised of a periodic pulse train; splitting the input signal into two optical signals; forming at least one cluster of pulses in each of the two optical signals, where amplitude and phase of pulses in each cluster varies amongst pulses; and coherently combining the two optical signals in the time domain using a resonant cavity and thereby forming an optical output signal having a solitary pulse for each cluster of pulses, where the two optical signals are combined in a manner such that pulse energy of the solitary pulse exceeds two times highest energy of any given pulse in the cluster of pulses from the two optical signals.
 2. The method of claim 1 wherein forming at least one cluster of pulses further comprises amplitude modulating the two optical signals and phase modulating the two optical signals.
 3. The method of claim 1 further comprises forming at least one cluster of pulses so that clusters of pulses in the two optical signals are complementary with each other in the spectral domain.
 4. The method of claim 1 further comprises coherently combining the two optical signals using a Fabry-Perot interferometer.
 5. A method for coherently combining pulses of a pulsed optical signal in the time domain, comprising: receiving an optical input signal comprised of a periodic pulse train; splitting the input signal into N optical signals; amplitude modulating each signal in the N optical signals to form clusters of pulses in each signal, such that amplitude of pulses in each cluster varies amongst pulses; phase modulating each signal in the N optical signals, such that phase of pulses in each cluster varies amongst pulses; and coherently combining each of the phase modulated signals in the time domain using one or more resonant cavities and thereby forming an optical output signal having a solitary pulse for each cluster of pulses, where the N optical signals are combined in a manner such that pulse energy of the solitary pulse exceeds N times highest energy of any given pulse in the cluster of pulses from the N optical signals.
 6. The method of claim 5 further comprises amplifying each of the phase modulated signals prior to coherently combining the phase modulated signals.
 7. The method of claim 5 further comprises phase locking each of the N optical signals.
 8. The method of claim 5 further comprises coherently combining each of the phase modulated signals using a plurality of Fabry-Perot interferometer cavities.
 9. The method of claim 8 further comprises arranging the plurality of Fabry-Perot interferometers in a parallel cascading manner.
 10. The method of claim 8 further comprises arranging the plurality of Fabry-Perot interferometers in a sequential cascading manner.
 11. A method for coherently combining pulses of a pulsed optical signal in the time domain, comprising: receiving an optical input signal comprised of a periodic pulse train; forming at least one cluster of pulses in the pulse train, where amplitude and phase of pulses in the at least one cluster varies amongst pulses; and coherently combining pulses in the at least one cluster in the time domain to form an optical signal with a solitary pulse using a resonant cavity, where the resonant cavity has an entirely reflective cavity.
 12. The method of claim 11 wherein forming at least one cluster of pulses further comprises amplitude modulating the optical input signal and phase modulating the optical input.
 13. The method of claim 11 further comprises coherently combining pulses using a Gires-Tournois interferometer cavity.
 14. A method for coherently combining pulses of a pulsed optical signal in the time domain, comprising: receiving an optical input signal comprised of a periodic pulse train; amplitude modulating the optical input signal to form at least one clusters of pulses, such that amplitude of pulses in each cluster varies amongst pulses; phase modulating the optical input signal, such that phase of pulses in each cluster varies amongst pulses; splitting the optical input signal into N optical signals; amplifying each of the N optical signals; coherently combining each of the amplified optical signals into a combined signal using a beam combiner; and coherently combining pulses in the at least one cluster of pulses in the combined signal in the time domain using a resonant cavity having an entirely reflective cavity, thereby forming an optical output signal having a solitary pulse for each cluster of pulses.
 15. The method of claim 14 further comprises phase locking each signal of the N optical signals.
 16. The method of claim 14 further comprises coherently combining pulses in each cluster of pulses using a Gires-Tournois interferometer.
 17. The method of claim 16 further comprises coherently combining pulses using N Gires-Tournois interferometers arranged in a sequential cascading manner, such that the energy increase of a solitary pulse is approximately 2N.
 18. The method of claim 16 further comprises coherently combining pulses using a plurality of Gires-Tournois interferometers arranged in two or more sequential stages, where each stage includes N Gires-Tournois interferometers arranged in a sequential cascading manner, where the Gires-Tournois interferometer in each subsequent cascade are of equal roundtrip length, which is at least 2N times longer than in the previous cascade. 